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Abstract 

Financial markets have developed a lot of strategies to control risks induced 
by market fluctuations. Mathematics has emerged as the leading discipline to 
address fundamental questions in finance as asset pricing model and hedging 
strategies. History began with the paradigm of zero-risk introduced by Black 
& Scholes stating that any random amount to be paid in the future may be 
replicated by a dynamical portfolio. In practice, the lack of information leads 
to ill-posed problems when model calibrating. The real world is more complex 
and new pricing and hedging methodologies have been necessary. This chal- 
lenging question has generated a deep and intensive academic research in the 
20 last years, based on super-replication (perfect or with respect to confidence 
level) and optimization. In the interplay between theory and practice, Monte 
Carlo methods have been revisited, new risk measures have been back-tested. 
These typical examples give some insights on how may be used mathematics 
in financial risk management. 
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1. Introduction 



Financial markets have become an important component of people's life. All 
sorts of media now provide us with a daily coverage on financial news from all mar- 
kets around the world. At the same time, not only large institutions but also more 
and more small private investors are taking an active part in financial trading. In 
particular, the e-business has led to an unprecedented increase in small investors 
direct trading. Given the magnitude of the potential impact a financial crisis can 
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have on the real side of the economy, large, and even small, breakdowns have re- 
ceived particularly focal attention. Needless to say, the rapid expansion of financial 
markets calls upon products and systems designed to help investors to manage their 
financial risks. The financial risk business now represents more than $15 trillion an- 
nually in notional. A large spectrum of simple contracts (futures, options, swaps, 
etc.) or more exotic financial products (credit derivatives, catastrophe bonds, ex- 
otic options, etc.) are offered to private investors who use them to transfer financial 
risks to specialized financial institutions in exchange for suitable compensation. A 
classic example is the call option, which provides a protection in case of a large 
increase in the underlying asset price . More generally, a derivative contract is an 
asset that delivers a payoff H(u>) at maturity date, depending upon the scenario io. 

As argued by Merton [TJJ, the development of the financial risk industry would 
not have been possible without the support of theoretical tools. Mathematics has 
emerged as the leading discipline to address fundamental questions in financial risk 
management, as asset pricing models and hedging strategies, based on daily (in- 
finitesimal) risk management. Mathematical finance, which relates to the applica- 
tion of the theory of probability and stochastic processes to finance, in particular the 
Brownian motion and martingale theory, stochastic control and partial differential 
equations, is now a field of research in its own right. 

2. The Black & Scholes paradigm of zero-risk 

It is surprising that the starting point of financial industry expansion is "the 
Brownian motion theory and Ito stochastic calculus" , first introduced in finance by 
Bachelier in this PhD thesis (1900, Paris), then used by Black, Scholes and Merton 
in 1973. Based on these advanced tools, they develop the totally new idea in the 
economic side that according to an optimal dynamic trading strategy, it is possible 
for option seller to deliver the contract at maturity without incurring any residual 
risk. At this stage, any people not familiar with stochastic analysis (as the majority 
of traders in the bank) may be discouraged. As Foellmer said in Bachelier Congress 
2000, it is possible to reduce technical difficulties, and to develop arguments which 
are essentially probability-free. 

A DYNAMICAL UNCERTAIN WORLD 

The uncertainty is modelled via a family £1 of scenarios uj, i.e. the possible 
trajectories of the asset prices in the future. Such paths are described as posi- 
tive continuous functions u> with coordinates X t — uj(t), such that the continuous 
quadratic variation exists: [A] t (cj) = lim„ X)t <t t eD ~ -^■u) 2 along the se- 

quence of dyadics partitions D n . The pathwise version of stochatic calculus yields 
to Ito's formula, 

f{t,X t )(u) = /(Mo) + Jo &{s,X.){w)dX.{u) + J* f t (a,X a ){u)dt 

1 A Call option provides its buyer with the right (and not the obligation) to purchase the risky 
asset at a pre-specified price (the exercise price) at or before a pre-specified date in the future 
(maturity date). The potential gain at maturity can therefore be written as (Xj- — K) + , where 
Xt denotes the value of the underlying asset at time T 
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+/,fj/»(a.)(«M. 

The second integral is well denned as a Lebesgue-Stieljes integral, while the first 
exists as Ito's integral, defined as limit of non-anticipating Ricmann sums, (where 
we put S t = F x (t,X t )), E ti <t,neDjtM(Xt i+1 -X ti ). 

From a financial point of view, the Ito's integral may be interpreted as the 
cumulative gain process of trading strategies: S t is the number of the shares held 
at time t then the increment in the Riemann sum is the price variation over the 
period. The non anticipating assumption corresponds to the financial requirement 
that the investment decisions are based only on the past prices observations. The 
residual wealth of the trader is invested only in cash, with yield rate (short rate) r t 
by time unit. The self-financing condition is expressed as: 

dV t = r t (V t - S t .X t )dt + S t .dX t = r t V t dt + S t . (dX t - r t X t dt), V = z. 

Hedging derivatives: a solvable target problem 

Let us come back to the problem of the trader having to pay at maturity T the 
amount in the scenario u — K) + for a Call option). This target 

has to be hedged (approached) in all scenarios by the wealth generated by a self- 
financing portfolio. The "miraculous" message is that, in Black & Scholes world, a 
perfect hedge is possible and easily computable, under the additional assumption: 
the short rate is deterministic and the quadratic variation is absolutely continuous 
d[X] t — a(t, X t ) Xfdt. The (regular) strictly positive function a(t,x) is a key 
parameter in financial markets, called the local volatility. 

Looking for the wealth as a function f(t,X t ), we see that, given Ito's formula 
and self-financing condition, 

df(t, x t ) = fl(t, x t )dt + f x (t, x t )dx t + y x ' x (t, X t )X?o- 2 (t, X t ) dt 

= f{t, X t )r t dt + S(t, Xt)(dX t - X t r t dt) 

By identifying the dX t terms (tanks to assumption a(t, x) > 0), S(t, X t ) = f x (t, X t ), 
and / should be solution of the following partial differential equation, Pricing PDE 
in short, 

fl(t,x) + y x ' x (t,x)x 2 a 2 (t,x) +f x (t,x)xr t - f(t,x)r t = 0, f(T,x) = h(x) (2.1) 

The derivative price at time to must be f(t ,x ), if not, it is easy to generated profit 
without bearing any risk (arbitrage). That is the rule of the unique price, which 
holds in a liquid market. 

The PDE's fundamental solution q(t,x,T,y) (h(x) = 5 y (x)) may be inter- 
preted in terms of Arrow-Debreu "states prices" density, introduced in 1953 by these 
Nobel Prize winners for a purely theoretical economical point of view and by com- 
pletely different arguments. The pricing rule becomes: f(t,x) = J h(y)q(t,x,T,y)dy. 
q is also called pricing kernel. When a(t, x) = at, the pricing kernel is the log-normal 
density, deduced from the Gaussian distribution by an explicit change of variable. 
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The closed formula for Call option price is the famous 2 Black and Scholes formula, 
which is known by any practitioner in finance. The impact of this methodology was 
so important that Black (who already died) , Scholes and Merton received the Nobel 
prize for economics in 1997. 

In 1995, B.Dupire [H] give a clever formulation for the dual PDE (one dimen- 
sional in state variable) satisfied by q(t,x,T,y) in the variables {T 7 y). If C(T,K) 
is the Call price with parameters (T,K) when market conditions are (to,Xo), then 



C' T {T,K) = ±<r 2 (T,K)K 2 C£ K (T,K)~rKC' K (T,K), C(t ,x ) = (x - K) 
In short, if r t = 0, C' T {T,K) = ^a 2 (T, K) K 2 C'^ K (T, K) 

3. Model calibration and Inverse problem 



In pratice, the main problem is the model calibration. Before discussing that, 
let me put the problem in a more classical framework. Following P. Levy, asset price 
dynamics may be represented through a Brownian motion via stochastic differential 
equation (SDE) 



where the Brownian motion W may be viewed as a standardized Gaussian noise 
with independent increments. 

The local expected return /i(t, X t ) is a trend parameter appearing for the first 
time in our propose. That is a key point in financial risk management. Since this 
parameter does not appear in the Pricing-PDE, the Call price does not depend 
on the market trend. It could seem surprising, since the first motivation of this 
financial product is to hedge the purchaser against underlying rises. By using 
dynamical hedging strategy, the trader (seller) may be also protected against this 
unfavorable evolution. For a statistical point of view, this point is very important, 
because this parameter is very difficult to estimate. 

In the B & S model with constant parameters, the volatility square is the 
variance by time unit of the log return ln(JT t ) — \n{X t _h) = R t . If the only avail- 
able information is given by asset price historical data, the statistical estimator to 
be used is the empirical variance, computed on a more or less long time period, 
a 2 = jrziYli^iRti ~ Rn) 2 , where Rn — jf^f^Rti- This estimator is called 
historical volatility. 

However, traders are reserved in using this estimator. Indeed, they argue that 
financial markets are not "statistically" stationary and that past is not enough to 

2 In the Black-Scholes model with constant coefficients, the Call option price C BS (t, x, K,T) is 
given via the Gaussian cumulative function M(z) = e 2~diy,and 9 = T — t, 



dX t - X t (ji{t, X t )dt + a(t, X t )dW t ), X t0 = x 




/A')) - Ke- r0 M(d o (e,x/K)) 
8, d 1 (8,x/K) = d (8,x/K) + a^8 



Moreover A(t,x) = d x C BS (t, x, t + 9, K) = JV"(di(0, x/K)). 
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explain the future. When it is possible, traders use additional information given by 
quoted option prices and translate it into volatility paramctrization. The implied 
volatility is defined as: C obs (T, K) = C BS (t ,x 0l T, K, a trn P). 

Moreover, the A of the replicating portfolio is A irn P = d x C BS (t , x ,T, K, a imp ). 

This strategy is used dynamically by defining the implied volatility and the 
associated A at any renegotiation dates. It was this specific use based on the 
hedging strategy that may explain the Black & Scholes formula success. 

This attractive methodology has been appeared limited: observed implied 
volatilities are depending on the option parameters (time to maturity and exer- 
cise price) (implied volatility surface) in complete contradiction with B&S model 
assumptions. In particular, the market quotes large movements (heavy tail dis- 
tribution) higher than in the log-normal framework. The first idea to take into 
account this empirical observation is to move to a model with local volatility 
a(t,x). The idea is especially attractive since the Dupire formula (2.2) gives a sim- 
ple relation between "quoted Call option prices" and local volatility: a 2 (K 7 T) = 
2C' T {T, K)(K 2 C'k K (T, K))~ l . The local volatility is computable from a continuum 
of coherent (without arbitrage) observed quoted prices. Unfortunately, the option 
market is typically limited to a relatively few different exercize prices and maturi- 
ties; a naive interpolation yields to irregularity and instability of the local volatility. 

Ill-posed inverse problem 

The problem of determining local volatility can be viewed as a function ap- 
proximation non linear problem from a finite data set. The data set is the value 
Cij of the solution at (to,£o) of Pricing-PDE with boundary conditions hij(x) — 
(x — Kij) + at maturity T^. Setting the problem as PDE's inverse problem yields to 
more robust calibration methods. These ideas appear for the first time in finance 
in 1997 but the problem is not classical because of the strongly non linearity 
between option prices and local volatility; the data set is related to a single given 
initial condition. 

Prices adjustment is made through a least square minimization program, in- 
cluding a penalization term related to the local volatility regularity. 

G(°) = Eij "id (/(*o, x ,hij,Ti, a(., .)) - C™ s )\ 
J (a, a) = a||Ver|| 2 + G{&) -> min CT . 

Existence and uniqueness of solution is only partially solved 0]. Using large devia- 
tion theory, the asymptotic in small time of local volatility is expressed in terms of 
implied volatility: a^ d (K, to)- 1 = ln(f /* ^f^. 

Avellaneda & alii [2] have used another penalization criterion based on a 
stochastic control approach ; the control is the volatility parameter itself constrained 
to be very close to a prior volatility (r)(a) = \a(t,x) — ao(t,x)\ 2 for instance). The 
gradient criterion is replaced by K(a) — U(to,XQ, a) where U(T, x, a) = and 

U' t {t,x) + \o- 2 {t,x) x 2 U'J. x {t 1 x) + rxU' x {t,x) - U(t,x) +r)(<r(t,x)) = 0. 



4. Portfolio, duality and incomplete market 
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In the previous framework, options market may be entirely explained by un- 
derlying prices. In economic theory it corresponds to market efficiency: a secu- 
rity price contains all the information on this particular security. In option world, 
the observed statistical memory of historical volatility leads naturally to consider 
stochastic volatility models with specific uncertainty 

dX t = Xt{ii(t,Xt,Yt)dt + *(t,X t ,Y f )dW}), dY t =r l {t,X u Y t )+ 1 (t,X u Y t )dW? 

where dW 1 and dW 2 are two correlated Brownian motions. 7 is the volatility of 
the volatility. What does it change ? In fact, everything ! Perfect replication by 
a portfolio is not possible any more ; the notion of unique price does not exist any 
longer... But, such a situation is often the general case. What kind of answer may 
wc bring to such a problem? 
Super-replication and Robust Hedging 

The option problem is still a target problem Or, to be replicated by a portfolio 

Vt(tt, S) = 7r + L Ei^^t depending on market assets X 1 . Constraints (size, 
sign...) may be imposed on investment decisions {8\) 3 . Let Vr be the set of all 
derivatives, replicable at time T by an admissible portfolio. Their price at to is the 
value of their replicating portfolio. 

Supcr-replicating Or is finding the smallest derivative Or G Vr which is greater 
than Ct in all scenarios. The super-replication price is the price of such a 
derivative. The C't replicating portfolio is the Ct robust hedging. 

There are several ways to characterise the super-replicating portfolio: 

1) Dynamic programming on level sets: this is the most direct (but least 
recent) approach. This method proposed by Soner & Touzi [21] has led the way for 
original works in geometry by giving a stochastic representation of a class of mean 
curvature type geometric equations. 

2) Duality: this second approach is based on the Vt "orthogonal space" , a set 
Qt of martingale measures to be characterised. The super-replication price is given 

by 

Co =sup QeQT Eg [Or]. 
Wc develop this last point, which is at the origin of many works. 
Martingale measures 

The idea of introducing a dual theory based on probability measures is due to 
Bachelier (1900), and above all to Harisson & Pliska (1987). The actual and achieved 
form is due to Delbaen & Schachermayer |2j and to the very active international 
group in Theoretical Mathematical Finance. 

A martingale measure is a probability such that: Wt G Vr, Eq[Vt] = Vq. Us- 
ing simple strategies (discrete times, randomly chosen), this property is equivalent 
to prices of fundamental assets (X[) are Q-(local) martingales: the best Xl +h - 
cstimatcd (w.r. to Q) given the past at time t is X\ itself. The financial game is 
fair with respect to martingale measures. 

3 For the sake of simplicity, interest rates are assumed to be null 
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When Vt contains all possible (path-dependent) derivatives, the market is said 
to be complete, and the set of martingale measures is reduced to a unique element 
Q, often called risk-neutral probability. This is the case in the previous framework. 
Dynamics become dX t = X t a(t, X t )dW® where W® is a Q-Brownian motion. This 
formalism is really efficient as it leads to the following path dependent derivative 
pricing rule: Cq — Eg (Or)- 

Computing the replicating portfolio is more complex. In the diffusion case, 
the price is a deterministic function of risk factors and the replicating portfolio only 
depends on partial derivatives. The general case will be mentioned in the paragraph 
dedicated to Monte-Carlo methods. 
Incomplete market 

The characterization of the set Qt is all the more delicate so since there are 
many different situations which may lead to market imperfections (non-tradable 
risks, transaction costs ...). 

An abstract theory of super-replication (and more generally of portfolio opti- 
mization under constraints) based on duality has been intensively developed. The 
super-replicating price process is showed |15 | to be a super-martingale with re- 
spect to any admissible martingale measure. Hence, by the generalization of the 
Doob-Meyer representation, the super-replicating portfolio is the 11 Qt- martingale" 
part of the super-price Kramkov-decomposition. 

Super-replication prices are often too expensive to be used in practice. How- 
ever, they give an upper bound to the set of possible prices. In the previously de- 
scribed stochastic volatility model, the super-replication price essentially depends 
on possible values of stochastic volatility: 

1. If the set is R + , then the super-replicating derivative of h(Xx) is h(Xx) where 
h is the concave envelop of h ; the replicating strategy is the trivial one: buying 
h' x (xo) stocks and holding them till maturity. 

2. If the volatility is bounded (up and down relatively to 0), the super-replication 
price is a (not depending on y) solution of 

h' t (t, x) + \ sup y (cr 2 (i, x, y)K x (t, x)) = 0, h(T, x) = h(x). 

When h is convex, h(t, x) is convex and the super-replication price is the one 
calculated with the upper volatility (in y). 

Calibration constraints may be easily taken into account without modifying this 
framework. We only have to assume that the terminal net cash flows of calibrating 
derivatives belong to Vt or equivalently we have to add linear constraints to the 
dual problem: C cal = sup{E Q (C T ) ; Q G Qt , Eg((X Ti - = C tJ }. 

Risk measures 

When super-replicating is too expensive, the trader has to measure his market 
risk exposure. The traditional measure is the variance of the replicating error. But 
a new criterion, taking into account extreme events, is now used, transforming the 
risk management at both quantitative and qualitative levels. 
Value at Risk 

The VaR criterion, corresponding to the maximal level of losses acceptable with 
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probability 95%, has taken a considerable importance for several years. Regulation 
Authorities have required a daily VaR computation of the global risky portfolio 
from financial institutions. Such a measure is important on the operational point 
of view, as it affects the provisions a bank has to hold to face market risks. VaR 
estimation (quantile estimation) and its links with extreme value theory |12j are 
widely debated in the market, just as by academics. 

Moreover, a huge debate has been introduced by academics P on the VaR 
efficiency as risk measure. For instance, its non-additive property enables banks to 
play with subsidiary creations. This debate has received an important echo from the 
professional world, which is possibly planning to review this risk measure criterion. 
Sub-additive and coherent risk measures are an average estimation of losses with 
respect to a probability family: p{X) — supq gQt Eq(— X). 

This characterisation has recently been extended to convex risk measures, by 
adding a penalization term depending on probability density (entropy for instance) . 

Risk measures and reserve price 

A trader willing to relax the super-replication assumption is naturally thinking 
in terms of potential losses with given probability (level confidence) of his replicat- 
ing error. It corresponds to the quantile hedging strategy. Other risk measures 
(quadratic, convex, entropic) may be used. Optimization theory is coming back 
when looking for the smallest portfolio, generating an acceptable loss. The ini- 
tial value of this portfolio is called the reserve price. Mean-variance and entropic 
problems have now a complete solution 13 . More surprisingly (because of non- 
convexity), this also holds for the quantile hedging problem AH these results 
are in fact sub-products of portfolio optimization in incomplete markets [01 or 

5. New research fields 

Monte-Carlo methods 

Dual version of super-replication problems, just as new risk measures, under- 
line the interest to compute very well and quickly quantities such as Eg (A) and 
more generally supg g g T Eq(X). For small dimensional diffusions, these quantities 
may be computed as the solution of some linear PDE (for the expected value) and 
non- linear PDE (for the sup). However, the computational efficiency falls rapidly 
with the dimension. That increases the interest for the so-called probabilistic meth- 
ods. 

The fundamental idea of Monte-Carlo methods is the computation of Eq(A) 
by simulation, i.e. by drawing a large number (N ~ 10 5 ) of independent scenarios 
u) 1 and taking the average value of the results -k J2i=i -^ito 1 )- Of course, this method 
does not work very well when being too naive, but convergence may be accelerated 
by different techniques. 

In the finance area, the important quantities are both the price and the sen- 
sitivities to different model parameters. Based on integration by parts, efficient 
methods have been developed to compute in a coherent manner prices and their 
derivatives |17|. In the case of path-dependent options, the derivative is taken with 
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respect to a Brownian motion perturbation (Malliavin Calculus). 

Very original is the actual research on solving, by Monte-Carlo methods, 
optimization problem expressed as "sup" of a family of expected values (super- 
replication prices). These solutions are based on dynamic programming, enabling 
to turn a maximization in expected value into a pathwise maximization. The for- 
mulation in terms of Backward SDE's, introduced by Peng and Pardoux in 1987 
and in finance 4 in we ^ describes this effect. 

In all cases, the problem is to compute conditional expectation by Monte 
Carlo, using the function approximation theory, or more generally random variable 
approximation in in the Wiener space (chaos decomposition). 

Problems related to the dimension, and statistical modelling 

Financial problems are usually multidimensional, but only few liquid financial 
products are depending on multi-assets. Even if the different market actors consider 
they can have a good knowledge of each individual asset behavior, the question is 
now to find a multidimensional distribution given each component distribution. This 
problem is a statistical one, known as the copula theory: Copula is a distribution 
function on [0, 1]™ with identity function as marginal. They are useful to give 
bounds to asset prices. Dynamically, the problem still to be solved is to find the 
local volatility matrix of multidimensional diffusion given the "Dupire" dynamics 
of each coordinate. 

High dimensional problems arise when computing bank portfolio VaR (with 
a number of observations less than that of risk factors), or hedging derivatives 
depending on a large number of underlying assets. Main risk factors may be very 
different in a Gaussian framework or heavy tail framework (Levy processes) 
Random matrix theory or other asymptotic tools may bring some new ideas to this 
question. 

By presenting the most important tools of the financial risk industry, I have 
voluntary left apart anything on financial asset statistical modelling, which may be 
the subject of a whole paper on its own. It is clear that the VaR criterion, the 
market imperfections are highly dependent on an accurate analysis of the real and 
historical world 0. Intense and very innovating research is now developed (High- 
frequency data, ARCH and GARCH processes, Levy processes with long memory, 
random cascades). 

6. Conclusion 

As a conclusion, applied mathematicians have been highly questioned by prob- 
lems coming from the financial risk industry. This is a very active world, rapidly 
evolving, in which theoretical thoughts have often immediate and practical fallout. 

4 A BSDE solution is a couple of adapted processes (Yt,Zt) such that 
-dY t = f(t, Y t ,Z t )dt - Z t dWt, Y T = C T . 
Thanks to comparison theorem, Yt = sup ig j Yj 1 is the solution of the BSDE with driver /(t, y, z) = 
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On the other hand, practical constraints raise new theoretical problems. This pa- 
per is far from being an exhaustive view of the financial problems. It is more a 
subjective view conditioned by my own experience. Many exciting problems, from 
both theoretical and practical points of view, have not been presented. May active 
researchers in these fields forgive me. 
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